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Abstract 

Charmless semileptonic decays, B — ► X u £u, are studied in a sample of 232 million BB decays 
recorded with the BaBar detector, in events where the decay of the second B meson is fully 
reconstructed. Inclusive charmless decays are selected in kinematic regions where the dominant 
background from semileptonic B decays to charm is reduced by requirements on the hadronic mass 
Mx and the momentum transfer q 2 . The partial branching fraction for B — ► X u lv decays for Mx < 
1.7 GeV/c 2 and q 2 > 8 GeV 2 /c 4 is measured to be AB(B -» X u £u) = (0.87 ± 0.09 stat ± 0.09 sys ± 
0.01th) x 10~ 3 . The CKM matrix element \V u b\ is determined by using theoretical calculations 
of phase space acceptances. Theoretical uncertainties in this extrapolation are reduced by using 
the inclusive b — > photon spectrum and moments of the b — > civ lepton energy and hadronic 
invariant mass. 
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1 Introduction 



The principal physics goal of the BaBar experiment is to establish CP violation in B mesons and to 
test whether the observed effects are consistent with the predictions of the Standard Model (SM). 
CP violating effects result in the SM from an irreducible phase in the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix which describes the couplings of the charged weak current to quarks. An improved 
determination of the magnitude of the matrix element |V^,|, the coupling strength of the b quark 
to the u quark, will contribute critically to tests of the consistency of the measured angles of the 
unitarity triangle of the CKM matrix. 

The precise determination of \V u b\ is very difficult as, due to the large charmed backgrounds, 
we are forced to measure partial branching fractions and extrapolate them to the full phase space 
by relying on QCD based theoretical calculations. These calculations make use of specific as- 
sumptions and are affected by different uncertainties. It is therefore important to make redundant 
measurements by using several experimental techniques, and different theoretical frameworks. Mea- 
surements of the inclusive and exclusive charmless semileptonic decays of B mesons are sensitive 
to different theoretical calculations, and therefore to different sources of systematic uncertainties. 
In addition, exploiting the available kinematic variables which discriminate between rare charmless 
semileptonic decays and the much more abundant decays involving charmed mesons, gives different 
sensitivities to the underlying theoretical calculations and assumptions. 

In inclusive measurements, three kinematic variables are discussed in the literature, each having 
its own advantages: the lepton energy (Eg), the hadronic invariant mass (Mj), and the leptonic 
invariant mass squared (q 2 ). The first measurements were restricted to the high end of the lepton 
spectrum where theoretical uncertainties are very large and therefore the extrapolation to the full 
spectrum becomes uncertain. Event selection based on Mx and q 2 allows us to select larger portions 
of phase space, but the underlying theoretical assumptions need to be carefully evaluated. 

Theoretical studies indicate that it is possible to reduce the theoretical error on the extrapolation 
by taking advantage of other kinematic variables or applying simultaneous cuts on Mx and q 2 in 
inclusive B — ► X u £u 4 decays [1]. In fact, while the Mx distribution has a large usable fraction of 
events, of the order of 70%, but depends on the shape function (SF) describing the Fermi motion of 
the b quark inside the B meson, the q 2 distribution is less sensitive to non-perturbative effects and 
less dependent on the calculation. Unfortunately, only a small fraction of events (about 20%) is 
usable with a pure q 2 selection. The study presented in [1] shows that a combined cut on Mx and 
q 2 may mitigate the drawbacks of the two methods while retaining good statistical and systematic 
sensitivities. 

BaBar has already published a determination of \V u b\ from a measurement of the inclusive 
charmless semileptonic branching fraction B(B — ► X u iv) [2] based on the study of the recoil of 
fully reconstructed B mesons and applying a kinematic cut on Mx< 1.55 GeV/c 2 , which resulted 
in the most precise determination of this quantity. Nevertheless the measurement was dominated 
(17% on the branching ratio, i.e. 8.5% on \V u b\) by the theoretical uncertainty on the underlying 
kinematic variable distributions, and therefore on the extrapolation to the full phase space. This 
result was obtained on a limited dataset, corresponding to about 80 fb _1 of integrated luminosity. 

Measurements of \V u b\ through inclusive charmless semileptonic decays B — > X u tv using a 
combination of Mx and q 2 ("Mx-<? 2 " analysis) have also been presented by BaBar at ICHEP 
2004 [3]. In this case, the extrapolation of the branching fraction measurements from a limited 
region of phase space to the full spectrum was done following Bauer et al [1] (BLL in the following). 

4 Charge- conjugate modes are implied throughout this paper. 
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The calculation for charmless semileptonic decays [4] implemented in the BaBar Monte Carlo 
simulation was also used to evaluate theoretical uncertainties. The preliminary result presented 
in this paper extends the analysis published in [3]. This new study is based on the same analysis 
strategy, but on a much larger dataset. 

This paper is organized as follows: Section 2 describes the detector, the data sample and the 
Monte Carlo simulation, including a description of the theoretical model on which our efficiency 
calculations are based. Section 3 describes the event reconstruction and selection, while in Sec- 
tion 4.1 the results of the Mx-q 2 analysis are presented. These results are interpreted by using the 
calculations from BLL and Neubert et al (BLNP) [5, 6, 7]. 

The SF parameters mb and /j, 2 described in Section 5 are an essential input for the determination 
of kinematic acceptances. In our previous result [3] we used two set of values for and /j, 2 , as 
determined from the analysis of the photon spectrum in B — > X s ~j decays from CLEO and Belle, 
respectively. For the present studies, we use the Belle result. The SF parameters can also be 
determined from the analysis of the electron energy and hadronic mass moments in b — ► civ decays. 
Therefore, in Section 5.1 we quote results based on the BaBar measurements of SF parameters 
with b — > civ moments. 

2 Data Sample and Simulation 

The data used in this paper were recorded with the BaBar detector [8] at the PEP-II collider. The 
total integrated luminosity of the data set is 210.7 fb _1 collected on the T(4S) resonance. The 
corresponding number of produced BB pairs is 232 million. We use Monte Carlo (MC) simulations 
of the BaBar detector based on GEANT [9] to optimize selection criteria and to determine signal 
efficiencies and background distributions. 

2.1 Simulation of B — > X u tv decays 

Charmless semileptonic B — > X u iv decays are simulated as a combination of both exclusive three- 
body decays to narrow resonances, X u = n, r/, p, uj, rj , and inclusive decays to non-resonant hadronic 
final states X u . The simulation of the inclusive charmless semileptonic B decays into hadronic states 
with masses larger than Im^ is based on a prescription by De Fazio and Neubert [4] (DFN), which 
calculates the triple differential decay rate, d 3 T / dq 2 dEg dsn (sr = M x ), up to 0(a s ) corrections. 
The motion of the b quark inside the B meson is incorporated in the DFN formalism by convolving 
the parton-level triple differential decay rate with a non-perturbative SF. The SF describes the 
distribution of the momentum k+ of the b quark inside the B meson. The two free parameters 
of the SF are A SF and \f F '. The first relates the B meson mass, mg, to the b quark mass, 
mf F = tub — A SF , and fi 2 = —Xf F is the average momentum squared of the b quark in the B 
meson. The SF parameterization used in the generator is of the form 

F(k+) = N(l-x) a e^ 1+a >, (1) 

where x = -j^p < 1 and a = —3(A SF ) 2 /\f F — 1. The first three moments of the SF must satisfy: 
A = 1, A x = and A 2 = -Af F /3. 

In the simulation the hadron system X u is produced with a non-resonant and continuous in- 
variant mass spectrum according to the DFN model. The fragmentation of the X u system into 
final state hadrons is performed by JETSET [10]. The exclusive charmless semileptonic decays are 
simulated using the ISGW2 model [11]. The resonant and non-resonant components are combined 
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such that the total branching fraction is consistent with the measured value [2] and that the inte- 
grated spectra agree with the prediction of Ref. [4] . In the evaluation of the associated uncertainty- 
all branching fractions and theory parameters are varied within their errors. 



2.2 Simulation of Background Processes 

To estimate the shape of the background distributions we make use of simulations of e + e~ — > 
T(45) — > BB with the B mesons decaying inclusively. The most relevant backgrounds are due 
to B — > X c £i> events. The simulation of these processes uses a Heavy Quark Effective Theory 
(HQET) parametrization of form factors for B — > D*£U [12], and models for B — ► D^-k(J> [13], 
and 5 — »■ D£u,D**il7 [11]. We also make use of a JETSET simulation of "continuum" e + e~ — »■ 
(q = u, d, s, c) events. 

3 Event Selection and Reconstruction 

The event selection and reconstruction and the measurements of branching fractions follow closely 
the strategy described in Ref. [2] and represent the basis for the M x -q 2 analysis presented here. 

In this paper we study the recoiled B candidate opposite of a fully reconstructed B in hadronic 
decay (B reco ), which is a moderately pure sample of B mesons with known flavor and four- 
momentum. We select -B reC o decays of the type B — > DY, where D refers to a charm meson, and Y 
represents a collection of hadrons with a total charge of ±1, composed of niir^+^K^+nsKg+n^ , 
where n\ + ri2 < 6, n% < 3, and < 3. Using D~ and D*~ (D° and D*°) as seeds for B° (B + ) 
decays, we reconstruct about 1000 different decay chains. Overall, we correctly reconstruct one 
B candidate in 0.3% (0.5%) of the B°B° (B + B~) events. The kinematic consistency of a _B reco 
candidate with a B meson decay is checked using two variables, the beam-energy-substituted mass 



fn.ES = \J s /^~Pb an d the energy difference, AE = Eb — \fs/^- Here -^/s refers to the total 
energy in the T(4S) center of mass frame, and ps and Eb denote the momentum and energy of the 
.B reco candidate in the same frame. For signal events the vies distribution peaks at the B meson 
mass, while AE is consistent with zero. We require AE = within approximately three standard 
deviations, as determined from each decay channel. 

A semileptonic decay of the other B meson (-B reC oii) is identified by the presence of a charged 
lepton with momentum in the -B reco ii rest frame (p*) above 1 GeV/c . In addition, the detection of 
missing energy and momentum in the event is taken as evidence for the presence of a neutrino. The 
hadronic system X is reconstructed from charged tracks and energy depositions in the calorimeter 
that are not associated with the -Breco candidate or the identified lepton. Care is taken to eliminate 
fake charged tracks, as well as low-energy beam-generated photons and energy depositions in the 
calorimeter from charged and neutral hadrons due to beam backgrounds. The neutrino four- 
momentum p v is estimated from the missing momentum four- vector p m iss = Pt(as) ~Pb icco ~Px ~Pe, 
where all momenta are measured in the laboratory frame, and Pt(4S) refers to the T(4S) meson. 

Undetected particles and measurement uncertainties affect the determination of the four- 
momenta of the X system and neutrino, and lead to a leakage of B — > X c lv background from 
the high M x into the low M x region (M x < 1.86 GeV/c 2 ). 

In the sample of reconstructed B decays (B rcco ) two backgrounds need to be considered: the 
combinatorial background from BB and continuum events, due to random association of tracks 
and neutral clusters, which does not peak in m-ES) and the BB background whose mES distribution 
has the same shape as the signal. After applying all selection criteria, the remaining combinatorial 
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background is subtracted by performing an unbinned likelihood fit to the m^s distribution. In 
this fit, the combinatorial background originating from e + e~ — ► qq (q = u, d, s, c) continuum and 
BB events is described by an empirical threshold function [14], and the signal is described by a 
modified Gaussian [15] peaked at the B meson mass. In addition, to further reduce the effects of 
the combinatorial background, only events with ttjes >5.27 GeV/c 2 are considered. 

To reject the background in the sample of semileptonic decays we require exactly one charged 
lepton with p* >lGeV/c, a total event charge of zero, and a missing mass consistent with zero 
( m miss < 0.5GeV 2 /c 4 ). These criteria partly suppress the dominant B — ► X c £u decays, many of 
which contain an additional neutrino or an undetected K® meson. 

In order to reject the peaking background coming from D* we use a dedicated method, the 
partial D* reconstruction. We explicitly veto the B° — > D* + l~V background by searching can- 
didates for such a decay with a partial reconstruction technique, that is only identifying the irf 
from the D* + — > D irf decay and the lepton: since the momentum of the 7r+ is almost collinear 
with the D* + momentum in the laboratory frame, we can approximate the energy of the D* + 
as E D *+ ~ m D »+ ■ E n J145 MeV/c 2 and estimate the neutrino mass as m 2 = (pb — Pd*+ ~ Pi) 2 - 
Events with m 2 > — 3 GeV 2 /c 4 are likely to be background events and are rejected. Finally, we 
veto events with charged or neutral kaons in the recoiling B to reduce the peaking background 
from B — > X c lv decays. Charged kaons are identified [8] with an efficiency varying between 60% 
at the highest (4 GeV/ c) and almost 100% at the lowest momenta. The pion misidentification rate 
is about 2%. The Kg — > 7t + tt~ decays are reconstructed with an efficiency of 80% from pairs of 
oppositely charged tracks with an invariant mass between 486 and 510 MeV/c 2 . 



4 Measurement of Charmless Semileptonic Branching Ratios 

To reduce the systematic uncertainties in the derivation of branching fractions, the observed number 
of signal events, corrected for peaking background and efficiency, is normalized to the total number 
of semileptonic decays B — > Xlv in the recoil of the -B reC o candidates. The number of observed 
-Breco events which contain a charged lepton with p* > 1 GeV/ c is denoted as N™ eas . It can be 
related to the true number of semileptonic decays, iV*[ Me and the remaining peaking background 
BG sU estimated with Monte Carlo simulation, by iV*[ ue = (N™ eas -BG sl )/efef = N sl /efef . Here 
ef refers to the efficiency for selecting a lepton from a semileptonic B decay with a momentum 
above p cu t in an event with a reconstructed B with efficiency ef. Figure 1 shows the result of the 
777.es fit used to determine N™ eas . 

If we denote as N™ eas the number of events fitted in the sample after all requirements, and 
with BG U the peaking background coming from semileptonic decays other than the signal, the true 
number of signal events N^ ue is related to them by 

M u - BG U - e sel e l e t M u , (2) 

where the signal efficiency e u sel accounts for all selection criteria applied on the sample after the 
requirement of a high momentum lepton. 

To measure BG U in the inclusive studies, the peaking background (BG U ) is estimated by per- 
forming a x 2 fit on the Mx-q 2 distributions, resulting from th-es fits in individual Mx-q 2 bins, with 
the shape of the background estimated from Monte Carlo simulation, and its normalization free to 
vary. 
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Figure 1: Fit to the ttt-es distribution for the sample with high momentum lepton. 



The ratio between the partial branching fractions for the signal and B — » Xlv decays is 
A£(signal) _ Nl rue _ {N~ - BG u )/{e u sel ) 
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B{B -» Xtv) 



jvrtrue 

Ajtrue 
si 
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-I e t 



meas 
si 



BG sl ) 



(3) 



The efficiency ratio 



is expected to be close to, but not equal to unity. Due to the difference in 



multiplicity and the different lepton momentum spectra, we expect the tag efficiency et and lepton 
efficiency e\ to be slightly different for the two classes of events, the largest effect coming from e/. 
The ratio was measured to be 1.204 db 0.033. The signal branching fraction is then obtained from 
R u / S i using the total semileptonic branching fraction oiB{B — > Xlv) = (10.83±0.19) %, which is the 
sum of the charm semileptonic branching ratio B(B — > X c lv) = (10.61 ±0.16(exp)±0.06(th))% [16] 
and the charmless semileptonic branching ratio B(B — > X u ti>) = (0.22±0.04(exp.) ±0.04(th))% [2], 
both measured in BaBar. 



4.1 Measurement of the Partial Branching Fraction 

Measurements done using only a Mx kinematic cut to reject the b — ► civ background are limited 
by the dependence on the SF. This can be overcome by selecting a phase space region where the 
SF effects are small, namely the region at large q 2 values [17]. In this way we find a trade-off 
between the statistical and theoretical uncertainties by loosening the Mx cut and applying a cut 
on q 2 . Moreover, since most of the theoretical uncertainties are due to the extrapolation from a 
selected kinematic region to the full phase space, measurements of partial branching fractions in 
different regions of phase space and their extrapolation to the full phase space can serve as tests of 
the theoretical calculations and models. 

In order to extract the partial charmless semileptonic branching ratio, AB(B — » X u £v), in a 
given region of the Mx-q 2 plane, we define as signal the events with true values of the kinematic 
variables in the chosen region, treating as background those that migrate from outside this region 
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Figure 2: Distributions of q 2 in two intervals of Mx- Points are data, the blue (medium shaded), 
magenta (dark shaded) and yellow (light shaded) histograms represent the fitted contributions 
from b — > ulv events inside true Mx< 1.7GeV/c 2 , q 2 > 8GeV 2 /c 4 , b — ► ulv events outside these 
requirements, and background events, respectively. 



because of the resolution. This means that in applying Eq. 3 we include the b — > ulv events outside 
the signal region in BG U and the quoted efficiencies refer only to events generated in the chosen 
(Mx-q 2 ) region. These efficiencies are computed on simulation based on the DFN model. However, 
the associated theoretical uncertainty on the final result is small compared to the extrapolation 
error to the full phase space. We divide the events into 32 non-equidistant two-dimensional bins of 
Mx and q 2 (4 bins in Mx and 8 in q 2 ), we fit the wes distribution to extract the yield in each bin, 
and we perform a two-dimensional binned fit of the entire Mx-q 2 distribution in order to extract 
the signal and background components. The result of the fit is shown in Fig. 2. We measure, out 
of 103590 ±474 background-subtracted semileptonic events (N™ eas — BG s i), 317 ±34 signal events 
(N™ eas - BG U ), above a background of 270 ± 5 events (BG U ). This, with e u sel = 0.319 ± 0.006, 
corresponds to a partial branching fraction in the signal region q 2 > 8GeV 2 /c 4 , Mx < 1.7GeV/c 2 
of: 

Ru/s\(B -»■ X u £v, M x < 1.7GeV/c 2 , q 2 > 8 GeV 2 /c 4 ) = 

= (0.80 ± 0.09 stat ± 0.09 syst ± 0.01 th ) x HT 2 , ( ' 

where the errors are due to statistics, experimental systematics and theoretical systematics, respec- 
tively. This gives the following value for the partial branching fraction: 

AB{B -> X u iv,M x < 1.7GeV/c 2 ,g 2 > 8GeV 2 /c 4 ) = , * 

= (0.87 ± 0.09 stat ± 0.09 syst ± 0.01 th ) x 10~ 3 . [ ' 

4.2 Systematic Uncertainties 

A breakdown of the systematic uncertainties is presented in Table 1. 

Uncertainties related to the reconstruction of charged tracks are determined by removing ran- 
domly a fraction of tracks corresponding to the uncertainty in the track finding efficiency (1.4% 
per track). 
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For photons, we correct for differences between data and MC in energy resolution, energy scale, 
and EMC crystal edge effects and assign the systematic uncertainty by repeating the analysis with- 
out applying the corrections. For single photon reconstruction no efficiency correction is applied, 
but a systematic uncertainty of 1.8% per photon is assigned. 

We estimate the systematic error due to particle identification (PID) by varying the electron 
and kaon identification efficiencies by ±2% and the muon identification efficiency by ±3%. The 
misidentification probabilities are varied by 15% for all particles. Effects due to Kl interactions 
have been estimated by removing all EMC deposits due to Kl when reconstructing Mj. 

The uncertainty of the B reco combinatorial background subtraction is estimated by changing 
the signal shape function to a Gaussian function instead of the empirical function of Ref. [15]. 
Furthermore, the parameters of the empirical function, which are kept fixed in the itt-es fits, are 
varied within their uncertainties. 

€ u € u 

The size of the Monte Carlo sample limits the accuracy on the determination of the ratio -^j-ti 
to 3%. 

The impact of the charm semileptonic branching fractions has been estimated by varying each 
of the exclusive branching fractions within one standard deviation of the current world average [18]. 
Similarly, the branching fractions of charm mesons for inclusive kaon production have been varied 
to estimate uncertainties in the kaon veto. 

To study the mixture of resonant decays among the charmless modes we also varied the number 
of charmless exclusive semileptonic decays by 30% for B — > ttIv and B — > p£u, by 40% for B — > 
loIu and by 100% for the remaining exclusive charmless semileptonic B decays. Using only the 
non-resonant model for the signal gives an estimate of the effects due to uncertainties in the 
hadronization model. Signal events where a gluon splits in a ss pair are varied by 30% in non- 
resonant events in order to obtain the associated systematic uncertainty. 

The uncertainties related to the knowledge of the SF are calculated by changing the functional 
form and varying the SF parameters as described in Section 5. 

5 Extraction of \V u b\ 

Using the Mx-q 2 analysis we measure the partial branching fraction for charmless semileptonic 
decays in a selected phase space region. To translate this into a measurement of the total branching 
fraction, and therefore |V^,|, we need the fraction of events inside the measurement region (referred 
to as "acceptance" in the rest of the paper) as an external input. 

In the following we use two different theoretical calculations of Bauer, Ligeti and Luke [17] 
(BLL) and Bosch, Lange, Neubert and Paz [5, 6, 7] (BLNP) for calculating acceptance corrections. 
Both BLL and BLNP use operator product expansions (OPE) to calculate QCD effects. 

5.1 Results using acceptances from BLL 

Bauer, Ligeti and Luke perform an OPE-based calculation to second order in the strong coupling 
constant a s and 6-quark mass rrif,. They focus on the region chosen for the measurement of the 
partial branching fraction where non-perturbative effects are small. In particular, they have shown 
that the theoretical uncertainties of the extrapolation to the full phase space are much reduced 
by restricting the selection to regions of higher values of q 2 , rather than just restricting Mx to a 
region below the charm meson mass. 
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Table 1: Systematic uncertainties in percent for the partial fraction AB(B — ► X u tv). 



Source 


cr(AB(B^X u £u)) 


Statistical error 


10.7 


Monte-Carlo statistics 


4.0 


Detector-related: 


Tracking efficiency 


2.2 


Neutral efficiency 


0.1 


Neutral corrections 


0.6 


K L 


2.0 


PID efficiency & misidentification 


2.5 


Detector uncertainties 


3.9 


-Brcco & fit: 


m E s fit 

e l e t 1 e l e t 


4.1 
3.0 


B SL branching ratios 


4.9 


D branching ratios 


0.1 


.Breco & fit errors 


7.1 


Signal: 


Composition of B — > decays 


3.0 


Hadronization 


3.0 


Gluon splitting to 


2.2 


Shape function parameters 


+1.6 
-1.9 


Shape function form 


0.3 


Signal efficiency 


5.1 


Total systematic error 


10.4 


Total error 


14.9 
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Based on these calculations we can convert the measured AB(B — > X u tv) into |V^j,| by 



_ 192vr3 ABjB^XJv) 
|Kfe '"V^^f G (6) 

where tb = 1.604 ± 0.012 ps [18] and G is a theoretical parameter calculated in the BLL ap- 
proach [17]. The first factor under the square root is 192-7T 3 / \TBG 2 F m\) = 0.00779. To extract 
\V-ub\, we take G = 0.27 as computed by BLL for mj,(lS) = 4.7 GeV/c 2 . We then infer the 6-quark 
mass in the IS scheme from the BaBar measurement of m\ m [16] by using the prescription in [19], 
obtaining mft(lS) = 4.74 GeV/c 2 . G is then recomputed by rescaling the original BLL value by 
the ratio (4.74/4.7) 9 [20], obtaining G = 0.291 ± 0.055. The 19% error on G, which turns into a 
9.5% error on \V u b\, is the sum in quadrature of uncertainties due to: residual SF effects, higher 
order terms in the a s perturbative expansion, a 80 MeV/c 2 uncertainty on the b quark mass, and 
C(ylQ CD /m 3 ) terms in the OPE expansion. The uncertainty on the b quark mass is the dominant 
source, contributing about 15% to the uncertainty on G. Eq. 6 yields 

\V ub \ = (4.82 ± 0.26 stat ± 0.25 syst ± 0.46 th+ s F ) x lO" 3 . (7) 
5.2 Results using the theoretical calculations by BLNP 

Bosch, Lange, Neubert, and Paz have performed calculations of the differential decay rates for 
B — ► X u lv and B — > X s j. The authors presented a systematic treatment of the SF effects, 
incorporating all known corrections to the rates, and provided an interpolation between regions 
of phase space that can be treated reliably by OPE calculations and others that depend on SF. 
They have introduced a parameterization of the SF. The parameters describing the SF cannot be 
calculated, rather they have to be taken from experiment. 

On the basis of these SF parameters, the partial rate for B — > X u lv can be predicted for the 
measured phase space, and related to |T4fc|> 

M^ ut 



Cut 



such that 



BLNP give results and uncertainties in terms of the reduced decay rate AC, defined in units of 
|K b | 2 ps- 1 . 

We rely on two measurements of these SF parameters, one based on the photon spectrum 
in B — > Xgj decays, the other on moments of the hadron mass and lepton energy spectrum in 
B — ► X c lv decays. The analysis of B — > X s ^y decays can be used to determine the SF parameters 
in a given renormalization scheme [21]. Likewise the moments of the lepton energy and hadronic 
invariant mass in B — > X c tv decays are sensitive to the heavy quark parameters, as shown in an 
OPE calculation [22] in the kinetic scheme. In both cases, the heavy quark parameters entering 
the calculations can be related to A SF and , see e.g. [23, 24]. 

The Belle Collaboration has measured the photon spectrum in B — > X s ~f decays [25] and, based 
on a fit to the spectrum, has determined A SF = 0.66GeV/c 2 and \f F = -0.40GeV 2 /c 4 [26]. They 
also provide a A\ 2 = 1 contour, which we use to estimate theoretical uncertainties. 
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These SF parameters translate to mf F = 4.52±0.07 and $ SF = 0.27±0.23 [27]. This results in 
A(" = (21.6±4.0±2'!)|T4b| 2 P s_1 > where the first error is due to the limited experimental knowledge 
of the SF parameters and the second to theory uncertainties, and consequently 

\V ub \ = (5.00 ± 0.27 stat ± 0.26 syst ± 0.46 SF ± 0.28th) x HT 3 , (10) 

where the errors are due to statistics, experimental systematics, shape function parameters and 
theoretical systematics, respectively. 

Alternatively, the BaBar collaboration has determined m£ m and fi 2,km in the kinetic mass 
scheme from fits to moments measured for B — > X c lv [16]. The values have been translated into 
the SF scheme by following the prescription in [23] resulting in mf F = 4.61 ± 0.08 GeV/c 2 and 

2 SF 

= 0.15 ± 0.07, with a correlation of -40%. The systematic error due to the uncertainty of the 
SF parameters is reduced, due to the significantly better precision obtained in the BaBar moments 
analysis. 

By using the results of the BaBar moments analysis we get A£ = (25.04 ±|'qq SF i 2. 45th) 
|^ub| 2 ps _1 . Again, the error is due to the limited experimental knowledge of the shape function 
parameters. This translates into 

\V ub \ = (4.65 ± 0.24 stat ± 0.24 syst +°;^ SF ± 0.23 th ) x 10~ 3 . (11) 
6 Conclusions 

We have presented a study of charmless semileptonic decays and a measurement of the \V ub \ CKM 
matrix element, by using the combined information of the Mx-q 2 distribution to discriminate signal 
and background and to minimize the theoretical uncertainties. We give a measurement of the 
partial branching fraction of charmless semileptonic decays AB(B — > XJLv) for Mx < 1.7GeV/c 2 
and q 2 > 8 GeV 2 /c 4 and, by taking kinematic acceptances from two theoretical calculations by BLL 
and BLNP, extract \V ub \. 

The measured partial branching fraction B{B — > X u tv) in the region limited by Mx < 
UGeV/c 2 , q 2 > 8GeV 2 /c 4 is 

AB(B -► XJu, M x < 1.7 GeV/c 2 , q 2 > 8 GeV 2 /c 4 ) = , , 

= (0.87 ±0.09 stat ±0.09 syst ± 0.01 th ) xl0~ 3 . ( > 

We extract the CKM matrix element \V ub \ using different approaches. With acceptances calculated 
using the BLL calculations, we obtain: 

|Kb| BLL = (4.82 ± 0.26 stat ± 0.25 syst ± 0.46 th+S F) x 10" 3 . (13) 

Using the partial decay models calculated in the BLNP approach and by taking the shape function 
parameters from the Belle photon spectrum in B — > X s ~f and the BaBar analysis of B — > X c lv 
moments, we find: 

|KbKV>x s7 = (5-00 ± 0.27 stat ± 0.26 syst ± 0.46 SF ± 0.28 th ) x lO" 3 , (14) 
\ V ^BaBar^cIu = (4-65 ±0.24 stat ±0.24 syst +^ SF ± 0.23 tll ) xlO- 3 , (15) 

where the errors are due to statistics, experimental systematics, SF and theoretical systematics, 
respectively. 
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In conclusion, the total error on \V u b\ is dominated by the experimental and theoretical uncer- 
tainties of the shape function. Our results of \V u b\ using the two different calculations of BLL and 
BLNP are consistent with each other. For the BLNP calculations the two sets of shape function 
parameters coming from a fit to the photon energy spectrum and to the B — ► X c lv moments are in 
good agreement and thus give consistent results on \V U \\. Results based on partial branching frac- 
tion from the lepton spectrum and q 2 using BLNP appear to be consistent with our measurement 
but somewhat lower [28]. 

7 ACKNOWLEDGMENTS 

We are grateful for the extraordinary contributions of our PEP-II colleagues in achieving the 
excellent luminosity and machine conditions that have made this work possible. The success of 
this project also relies critically on the expertise and dedication of the computing organizations 
that support BaBar. The collaborating institutions wish to thank SLAC for its support and the 
kind hospitality extended to them. This work is supported by the US Department of Energy and 
National Science Foundation, the Natural Sciences and Engineering Research Council (Canada), 
Institute of High Energy Physics (China), the Commissariat a l'Energie Atomique and Institut 
National de Physique Nucleaire et de Physique des Particules (France), the Bundesministerium fur 
Bildung und Forschung and Deutsche Forschungsgemeinschaft (Germany), the Istituto Nazionalc 
di Fisica Nucleare (Italy) , the Foundation for Fundamental Research on Matter (The Netherlands) , 
the Research Council of Norway, the Ministry of Science and Technology of the Russian Federation, 
and the Particle Physics and Astronomy Research Council (United Kingdom). Individuals have 
received support from CONACyT (Mexico), the A. P. Sloan Foundation, the Research Corporation, 
and the Alexander von Humboldt Foundation. Finally, we would like to thank the many theorists 
with whom we have had valuable discussions, and further thank M. Neubert, B. Lange and G.Paz 
for making available for our use a computer code implementing their calculations. 

References 

[1] C. W. Bauer, Z. Ligeti and M. E. Luke, Phys. Rev. D64, 113004 (2001), [hep-ph/0107074]. 

[2] BaBar Collaboration, B. Aubert et al, Phys. Rev. Lett. 92, 071802 (2004), [hep-ex/0307062]. 

[3] BaBar Collaboration, B. Aubert et al, "Study of b — > ulv Decays on the Recoil of Fully 
Reconstructed B Mesons and Determination of |V^,|", [hep-ex/0408068]. 

[4] F. De Fazio and M. Neubert, JHEP 06, 017 (1999), [hep-ph/9905351]. 

[5] B. O. Lange, M. Neubert and G. Paz, "Theory of Charmless Inclusive B Decays and the 
Extraction of \V u b\\ [hep-ph/0504071]. 

[6] S. W. Bosch, B. O. Lange, M. Neubert and G. Paz, Nucl. Phys. B699, 335 (2004), [hep- 
ph/0402094]. 

[7] S. W. Bosch, M. Neubert and G. Paz, JHEP 11, 073 (2004), [hep-ph/0409115]. 

[8] BaBar Collaboration, B. Aubert et al, Nucl. Instrum. Meth. A479, 1 (2002), [hep- 
ex/0105044]. 



18 



[9] GEANT4 Collaboration, S. Agostinelli et ai, Nucl. Instrum. Meth. A506, 250 (2003). 

[10] T. SjSstrand, Comput. Phys. Commun. 82, 74 (1994). 

[11] D. Scora and N. Isgur, Phys. Rev. D52, 2783 (1995), [hep-ph/9503486]. 

[12] CLEO Collaboration, J. E. Duboscq et ai, Phys. Rev. Lett. 76, 3898 (1996). 

[13] J. L. Goity and W. Roberts, Phys. Rev. D51, 3459 (1995), [hep-ph/9406236]. 

[14] ARGUS Collaboration, H. Albrecht et ai, Phys. Lett. B318, 397 (1993). 

[15] Crystal Ball Collaboration, T. Skwarnicki, "A Study of the Radiative Cascade Transitions 
Between the V and T Resonances", DESY F31-86-02, Ph.D. thesis. 

[16] BaBar Collaboration, B. Aubert et ai, Phys. Rev. Lett 93, 011803 (2004), [hep-ex/0404017]. 

[17] C. W. Bauer, Z. Ligeti and M. Luke, "Precision Determination of |F u6 |", [hep-ph/0111387]. 

[18] Particle Data Group, S. Eidelman et ai, Phys. Lett. B592, 1 (2004). 

[19] M. Battaglia et ai, "The CKM Matrix and the Unitarity Triangle", [hep-ph/0304132]. 

[20] Z. Ligeti, private communication. 

[21] A. L. Kagan and M. Neubert, Eur. Phys. J. C7, 5 (1999), [hep-ph/9805303] . 

[22] P. Gambino and N. Uraltsev, Eur. Phys. J. C34, 181 (2004), [hep-ph/0401063]. 

[23] M. Neubert, Phys. Lett. B612, 13 (2005), [hep-ph/0412241], and private communications. 

[24] D. Benson, I. I. Bigi and N. Uraltsev, Nucl. Phys. B710, 371 (2005), [hep-ph/0410080]. 

[25] Belle Collaboration, P. Koppenburg et ai, Phys. Rev. Lett. 93, 061803 (2004), [hep- 
ex/0403004]. 

[26] A. Limosani and T. Nozaki, "Extraction of the 6-quark Shape Function Parameters Using the 
Belle B -> X s ~f Photon Energy Spectrum", [hep-ex/0407052]. 

[27] Belle Collaboration, I. Bizjak et ai, "Measurement of the Inclusive Charmless Semileptonic 
Partial Branching Fraction of B Mesons and Determination of | V u b | Using the Full Reconstruc- 
tion Tag", [hep-ex/0505088]. 

[28] BaBar Collaboration, B. Aubert, "Determination of \V u b\ from Measurements of the Electron 
and Neutrino Momenta in Inclusive Semileptonic B Decays", [hep-ex/0506036], submitted to 
Physical Review Letters. 



19 



